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Abstract 

We study the — Vs oscillation effects in the near detector of the MINOS experi- 
ment. Conceptually, the MINOS search for sterile neutrinos with mass ~ 1 eV realizes 
an interesting situation of partial decoherence of the neutrino state at the production. 
This corresponds to a difference of energies of the two mass eigenstates that is com- 
parable with or bigger than the width of the initial state (pion) . We show that these 
effects modify the MINOS bound on mixing of sterile neutrino for Amli ^ 0.5 eV^ 
and make the experiment insensitive to oscillations with Aml^ > 15 eV^. Oscillations 
with Amli = (1 — 3) eV^ could explain some deficit of events observed in the low 
energy bins in the near detector and correspondingly the excess of events in the far 
detector. 



1 Introduction 

Recently the MINOS collaboration has published a very stringent limit on the active to 
sterile neutrino mixing from neutral-current interactions pQ. The z/^ — z/^ mixing angle, 624, 
should be smaller than 

024 < 7° (90% C.L.) for > mg (1) 

at 1-3 mixing value, ^13 = 0. For maximally allowed 613 the limit becomes slightly weaker: 
024 < 8°. Eq. ([1]) corresponds to |t//^4p < 0.015. Being valid, as we will confirm, for 
Am^ < 0.5 eV^ this limit essentially excludes the sterile neutrino interpretation of the 
LSND [2j| and MiniBooNE [3] results. Indeed, the required mixing for their explanation 
increases when Am^ decreases, whereas the bounds on |f/;x4p from MINOS and |f/e4p from 
reactor experiments change weakly with Am^ in the range of interest. 

The limit in Eq. ([T]), however, can not be applied to oscillations with Am^ ~ 1 eV^ 
relevant for LSND and MiniBooNE. Indeed, the analysis in [T] has been performed assuming 
that no oscillation-induced change of the neutrino event rate is measurable in the near 
detector (ND). The ND is located at the distance L = 1.04 km from the beam target, and 
therefore the oscillations with Am^ ~ 1 eV^ can not be neglected. If for instance, the 
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neutrino energy equals = 1 GeV, the oscillation length is 2.5 km. The first oscillation 
minimum is at = 0.8 GeV0. 

Admittedly, the systematic errors in the ND are bigger than the oscillation effect and, 
therefore, the latter is not "measurable" by this detector alone. However, the search of 
oscillations and the MINOS bound on the angle 624 are based on the comparison of the 
signals in the ND and far detector (FD). Namely, the energy spectrum of the neutral current 
(NC) events in the ND is extrapolated to the FD and confronted with the data. This allows 
one to substantially reduce the systematic errors. It is the difference of the oscillation 
effects in the ND and FD that allows one to obtain the bound on the oscillation parameters. 
Therefore, oscillation effects in the ND can not be neglected when "propagating" oscillation 
predictions to the FD. 

The search for neutrino oscillations from free pion decay with /\m\2 ~ 1 eV^ realizes the 
conceptually interesting situation in which the coherence of the neutrino state is partially 
broken at the production. Recall that in a given experimental setup, coherence is destroyed 
at production if it is possible to identify in principle (using kinematics of the process) which 
mass eigenstate is produced (see |7j and refernces therein). 

For free pion decay and undetected muon the energy uncertainty, cr^;, is determined by 
the pion decay rate aE ~ ^-k- The difference of energies of the mass states equals ~ Am? /2E 
and therefore the parameter 

Aw? ^ Idee /^x 

f = = 27r— (2) 

can be considered as a measure of the decoherence of the neutrino state. We will call it 
the decoherence parameter. Here l^ee = l/F is the decay length and ly = At^E/ Am? is the 
oscillation length. If ^ ^ 1 the two mass states can be resolved whereas for ^ <C 1 the 
uncertanty is large and decoherence can be neglected. In the case of MINOS, due to the 
large size of the pipe-line {Ip = 675 m), pions undergo free decay if collisions are neglected. 
For instance, for E.„ ~ 10 GeV the decay length, l^ec = ctq-^t, ~ 560 m, is smaller than Ip. 
If Am^ ~ 1 eV^, one has ^ ~ 1 and decoherence at the production can not be neglected. 
This is the case for most of the pion spectrum in MINOS, namely, the size of the region 
of coherent neutrino production is given by the decay length and is comparable with the 
oscillation length. For larger energies of the pion the decay region is determined by the size 
of pipe. 

In this paper we show that partial decoherence of the neutrino state at the production 
modifies the standard oscillation result in the MINOS experiment. We find the spectra of 
events in the ND and in the FD with t/^ — Vs oscillations taken into account. We estimate 
the modification of the bounds on sterile neutrino mixing when oscillations at the ND are 
included. We notice that the deficit of the events in the low energy bins in the ND and 
some excess of events in the FD could be explained by oscillations into sterile neutrinos. 

""^The possibility of a small oscillation effect in the ND in MINOS for Am^ ^ leV^ was mentioned 
previously in [4] 
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We elaborate on the coherence condition in other experiments searching for sterile neutri- 
nos with masses around 1 eV (MiniBooNE, new proposed experiments) to which similar 

considerations can be applied. 

The paper is organized as follows. In Sec. 2 we consider the oscillation effects in the 
MINOS ND assuming that neutrinos are produced incoherently. In Sec. 3 we discuss the loss 
of coherence at the production and compute the effective oscillation probability considering 
coherent neutrino production along the pion trajectory. We will show that in the MINOS 
setup, the same result as in the incoherent case can be obtained from the wave packet 
consideration for free pion decay under some approximations. In Sec. 4 we compute the 
numbers of events in the MINOS ND and FD and discuss modifications of the bounds when 
oscillations in the ND are taken into account. Conclusions are presented in Sec. 5 



2 Oscillation effects in the near detector. The case of 
incoherent production 

In this section we compute the oscillation effect in the ND considering incoherent production 
of neutrinos along the pion trajectory. We find first the oscillation probability for neutrinos 
produced in certain space-time point {xs, ts) of the production region (decay pipe) and then 
integrate this probability folded with the number of pions decaying in the point {xs, ts) over 
the decay pipe. 

Let us consider pions with energy E^^ produced in the target at origin, a; = 0, and 
compute the neutrino flux obtained at a detector located at x = L taking into account 
oscillations. The flux of neutrinos is given by 

F,{E,)= f dxP^^{E,,L-x)F,{E,,x), (3) 

where P^^ is the — survival probability and Ip is the length of the decay pipe. 
Fy{Ey^ x)dx is the density of neutrino flux in the interval {x, x + dx), produced from pion de- 
cay at a distance x from the target. After appropriate integration over the angular variables, 
the neutrino flux is equal to 

poo 

F4E,,x)= / dE^F4E^)e-^^^^^''T{E^)K{E^,E,), (4) 

TV 

where Ft^[Et^) is the flux of pions with energy Et^^V — mT^Fo/ET^ is the decay rate in the 
laboratory frame, Tq — I/tq — 2.5 • 10~^ eV is the decay rate in the pion rest frame; 

ET ^ ^ (5) 
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is the minimal pion energy required to produce a neutrino with energy Ei,. Eq is the neutrino 
energy in the rest frame of pion 

2 2 

mt — m,, 

where we have neglected the mass of the neutrino. K{E.,^, Ey) is the probability that a pion 
with energy Et^ emits a neutrino with energy E^. Similarly, one could write the contribution 
to the neutrino flux from two-body K decay just by substituting m^^ — )■ mx-, — t- T^- For 
illustration purposes we will only consider here neutrinos from pion decay which, in fact, 
constitute by far the main component of the neutrino beam. 

Notice that the integral in Eq. (jlj) is nothing but the integral of the oscillation proba- 
bility multiplied by the density of pions over the production region, that is the incoherent 
summation over neutrino sources. This is realized when the size of neutrino wavepackets is 
much smaller than the typical size of the source, ax ^ Ip- 

In what follows for simplicity we assume that the energy of neutrinos detected in the 
ND is uniquely related to the pion energy: 

= aE^, (7) 

where a = const. This approximation allows us to elucidate the physics involved and still 
obtain rather precise results. We consider validity of the approximation and determine value 
of the parameter a in the Appendix. In the approximation of Eq. ([7j) 

K{E^,E,) = 6{E^-a-'E,) (8) 

and the integration in Eq. (jl]) is trivial giving 

E,{E,,x) = F^(a-iE,)e-r(^-)^r(^,), (9) 

where 

m,) = (10) 

Integration of F^{E^,x) over x, which corresponds to setting P^^ = 1 in Eq. (E]), gives the 
neutrino flux at the ND without oscillations: 

= F^{a-'E,) [1 - e-r(i^.)'p] . (n) 

This relation determines the flux of pions, F^, in terms of the non-oscillating neutrino flux. 
Inserting this pion flux into Eq. ([9]) we obtain 

V(E ]p-'^(E'')^ 

ME.,x) = F'UE.) • (12) 
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Then, according to Eq. ([3]), the z/^— flux in the presence of oscillations equal 

^'^(^-) = 1 t dxTiE,)e-^^^'^^^P,,iE,, L-x). (13) 

This expression has been obtained essentially in the factorization approximation: the oscil- 
lation probability is multiplied by the production probability and later it is multiplied by 
the detection probability (cross-section). 

The oscillation probability can be written as 

P^^{E,, L - x) = P + ^ sin^ 2^24 cos[er(L - x)], (14) 

where 

P = 1 - - sin^ 2^24 (15) 

is the averaged oscillation probability, and the decoherence parameter ^ introduced in ([2]) 
can be rewritten as 

Notice that ^ does not depend on the baseline or neutrino energy but only on intrinsic 
parameters of the pion and the neutrinos. 

Using the flux in Eq. ( 1T3|) we can introduce the ratio of fluxes with and without oscilla- 
tions which can be considered as the effective survival probability. From Eqs. (|T3|) and ([Hj) 
we obtain 

P,;^ = ^ = P + 1 sin^ 2^24^3^ ^ ' dxe-"^"" cos[er(L - x)], (17) 



and the integration in Eq. (]T71) gives 
sin^ 2,9 1 

Peff = P+ ^^^^ Jf^ ■ ^ _ ^_r;^ [cos 0L + ^ slu 0L - C""'' (cOs(0L " 0p) + ^ sin(0L - 0p 



(18) 

Here 

^ = '^^ ^ = ^^^^ ^^^^ 

are the oscillation phases acquired over the baseline L and the pipe line, Ip. 

The most noticeable deviation of Eq. f ITS]) from the standard formula for oscillations 
with the baseline L, 

Pstand = P+\ sin' 2^24 COS (/-L, (20) 

is related to the factors which depend on ^ that lead to a suppression of the depth of oscilla- 
tions and to a modiflcation of the oscillatory term. Eq. ( 1T8|) provides a clear interpretation 
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of the decoherence parameter. For very short decay pipe, Tip <^ 1, the standard oscillation 
formula with baseline L, Eq. f l20p . which corresponds to a fully coherent superposition of the 
mass eigenstates, is recovered in the limit ^ -C 1. On the other hand, in the limit ^ — )■ cxd, 
Peff reduces to the averaged probability which reflects complete decoherence. Therefore, in- 
deed ^ is a measure of the decoherence introduced in the production process of the neutrino. 
The decoherence effect thus disappears at large Fq or small Am^. Other useful expression 
for C, is 

In the limit of low energies, IpT ^ 1, we obtain 

sin^ 26*24 . , , t ■ a \ /oo\ 
F^"^ 2(TTe) + ^sm0i), (22) 

and the oscillation pattern does not disappear in spite of the integration over large distances. 
The reason is that with decrease of energy the oscillation length decreases as Zj, oc E^, and 
the effective region of the pion decay, F oc -E,^ decrease in the same way. 

In Fig. [T] we show the effective probability Peff, Eq. f lTS]) . as a function of neutrino 
energy for different values of Am|2- For comparison we also show the standard probability 
with baseline L. A few comments are in order. 



1. Due to the large size of the neutrino production region, (and the effect of partial 
decoherence at the production) positions of the first and higher oscillation minima are 
shifted to lower energies in comparison with pointlike production case with baseline 
L. The relative shift increases with Am^g and decreases with energy. 

2. The depth of oscillations is suppressed due to the decoherence effect in comparison 
to the standard oscillation case. This suppression becomes stronger with increase of 
Am|2 and therefore ^, and decrease of energy. 

3. For characteristic neutrino energy = 2 GeV (responsible for the observed events 
with energy ~ 1 GeV) and Aml2 = 1 eV^ the oscillation effect in the ND is Pnd — 
0.97, as compared to the average suppression, Ppjj = 0.95, in the FD. Therefore 
Pfd/Pnd = 0.98 instead of 0.95 without oscillations in ND. So, clearly for this Aml2 
the bound on mixing parameter should be substantially weaker. For Aml2 = 0.5 eV^ 
the oscillation effect is 0.008 and the corresponding ratio, Pfd/Pnd = 0.96, is close 
to the result which MINOS uses. 



6 



F/Fo F/Fo 




Figure 1: The ratio of neutrino fluxes with and without oscillations as a function of neutrino 
energy for different values of Am\2'- 0.5, 1 , 2, 4 eV^, and sin^ 2^24 = 0.1. For comparison 
the standard oscillation probability with baseline L = 1.04 km is shown with dashed lines. 



3 Coherence loss at the production 

In this section, we find the oscillation effect assuming that neutrinos are coherently emitted 
along the pion trajectory. We will integrate the amplitude of oscillation over the production 
region and then compute the probability of the process. 

But before that, let us elaborate more on the decoherence parameter ^. According to 
the exact energy-momentum conservation in the decay process the difference of energies of 
the neutrino mass eigenstates equals in the pion rest frame 

- ^- (23) 
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This difference is comparable with the width 

AE% ~ To (24) 
at Am42 ~ 7 eV^. The ratio of the two parameters, 



at (25) 



. n n.m2n "I'n "T '"'II , X 

2 = ■ f • 26 



To 2m^ro 

coincides with the coherence parameter up to factor a ~ 0.45. Notice that the difference of 
momenta of the eigenstates is even bigger: 

2m^ m\ — 
Numerically, 

For ^ ^ 1 the energy splitting is larger than the uncertainty in energy, and therefore the 
neutrino mass eigenstates are produced incoherently. 

This consideration practically does not depend on the reference frame since violation 
of coherence is a physical phenomenon. In a frame where the pion moves with the Lorentz 
factor 7^ and neutrinos are emitted in the direction of pion we have Fq — )■ 27^Fo, and 

2 

AE°2 ^ ll^AEl^ 2 " 2 ' (28) 
mt — mi 

TV fl 

SO that the ratio changes roughly by a factor of 2 

If ^ ^ 1 the mass eigenstates are produced in the coherent state. In the MINOS case we 
deal with the intermediate situation when results depend on the shape of the wave packet 
of neutrino. Indeed, the expression in Eq. flTS]) depends on the exponential decay factor, 
which turns out to be the shape factor of the neutrino wave packet (squared) from free pion 
decay. 

In a reference frame where the pion has energy the decay length equals Idee = 
£',r/(^7rFo) = E^/larrLTrTo), and the oscillation length of neutrinos = AttE^/ Am\2. The 
ratio of the two lengths 

Ifjer Am^ f , , 

Ji££ = = _2_ ^29) 

li, 47ram^Fo 27r 

equals up to the factor l/27r the decoherence parameter. Therefore the condition of coher- 
ence is Idee ^ 

Consistent treatment of decoherence effects can be done in terms of wave packets. The 
produced z/^ state can be described as 



and for the detected state we have 



\v,{L - x)) =Y.U;,g^{x - L)e^^l(^-^)|i..). (31) 

k 

Here, \l/j(x, t) is the wave function of the mass eigenstate and z/j, gf and are the shape 
factors of the neutrino wave packets which correspond to the production and detection of 
the ith mass eigenstate. Also, pi and p[ are the mean momenta of the these wave packets, 
Ei = \/pf + mf and Vi is the group velocity of ith mass eigenstate (see [7] for details). 

We now show that under certain assumptions the wave packet calculation leads to the 
same result as in Sec. 2. In the case of neutrinos produced in pion decay the neutrino wave 
function in the space-time point {x,t), as well as the shape factor g^, can be found in the 
following way. The neutrino wave packet is formed as a result of an integration over the 
production region of partial plane waves emitted from each point of the region: 

<i/i{x, t) = j dxs j c/t5e-^^*«e'P'("-"«)-'^'(*-*s). (32) 

Here e~5'"*s gives the amplitude of probability to find the pion at the moment of time ts- 
Due to strong interactions the pion is produced in a small spacetime region around a; = 0, 
t = with very short wave packets. The production region is thus approximately described 
by a delta function 

5{xs-VT,ts)- (33) 

Furthermore, because the baseline is small for the ND, the spread of the wave packets due 
to the difference in group velocities is negligible. Therefore, the shape of the neutrino wave 
packet (exponential in this case) is conserved. The conservation of the shape implies in turn 
that there is a one-to-one correspondence between the points of the wave packet in a given 
point {x,t) and in point {xs.ts) at the production. This correspondence is expressed as 

Xs x-xs 

t ~ \ . (34) 

Here we take into account that before neutrino appears in the point x, t pion travels distance 
Xs from the origin with velocity w^, and then neutrino travels distance x — xs with average 
group velocity Vy. 

From Eqs. and (l33l) we obtain 

v^jv^t - x) , . ,x Xs jvut - x) 
xs{x,t) = , ts{x,t) = — = , (35) 

and the second equality leads to another delta function for ts'- 

s(t,-^^). (36) 
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The two deltas, Eqs. ( 15^ and (15^ . remove the integration over xs and ts in (152]) giving 



where and are now functions of x and t as in Eq. fl35|l . Therefore, according to Eq. fl30|) . 
the shape factor equals 

c/f (x - = e^^^*^e-^P^^«+^-^'*s^ ts = ts{x,t) oc - x). (38) 

The amplitude of probability of the z/^ detection is given by 

^mmW = j dxML~x)Mx,t)) (39) 
= ^\U^i\'' I dxgf{x-v,t)gf\x-L)e'P^''-'''^'e'^P^-P'^^^^^^^ (40) 

Since the size of the ND is much smaller than the decay pipe (production region) we can 
take g^{x — L) oc 6{x — L) and therefore 

A,Ut) = E \U,^\'9fiL - v.t)e^^^'^-^^''. (41) 

i 

The probability equals 

/ + 00 P + OO 

dt\A^,\^= / dtE(?f(L-t;,t)(7f (L-t;fct)|f/^,nf/^fc|V<^-(^'*), (42) 
i,k 

where (f)^f^{L,t) = {pi — Pk)L — {Ei — Ek)t = ApL — AEt . Using the expressions for the 
shape factors in Eq. (1551) we obtain 

/ + 00 
dte-^ts{L,t) J2 |f/^,nf/^,|2e*'/'»^(^'*), (43) 
■~ i,k 

where 

(l)ikiL,t)=Ap[L-v^tsiL,t)] - AE[t-tsiL,t)] (44) 

is the total phase difference between the mass eigenstates arriving at the detector in the 
moment of time t. The time factor in the last term of Eq. (H4|) can be rewritten as 



Then using the relation 



AE ^ Vi^Ap + ^Am^ (45) 
2E 
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(which follows from the dispersion relation) we obtain 



(L - vJs) 



(46) 



2E 



which exactly coincides with the standard oscillation phase. Inserting this expression into 
Eq. (1431) and changing the integration variable t to x = xs = v-jrts, which varies in the limits 
— Ip, we find 



This coincides (after normalization) with the expression found in Eq. (1131) . Thus, the result 
for completely coherent neutrino emission along the whole neutrino trajectory coincides with 
the case of completely incoherent emission which corresponds to very short wave packet of 
neutrino. Essentially, here we have a situation that resembles a coordinate-space version of 
the Kiers-Nussinov- Weiss theorem [8] where a fully incoherent ensemble of neutrino states 
is physically indistinguishable from neutrinos produced coherently if both cases lead to the 
same energy distribution. 

If pions collide, the shape factor of the neutrino wave packet is not determined by the 
decay function of pion [9J. The coherent emission will occur only between two consecutive 
collisions with substantial momentum transfer. This usually leads to the Gaussian shape 
factor for the neutrino wave packet with the width determined by mean free path of the 
pion. The consideration is simple if ax is much smaller than the decay length: cr^. ^ I decay 
Suppose the pion emits a neutrino coherently in some region of size centered at xs- 
Then, instead of leading to Eq. (1141) . the integration over the emission region gives for the 
oscillation probability 



where Siodo'x/li/) is the localization factor. With collisions, it is this probability that should 
be used in computations of the neutrino flux in Eq. (|T3|) and the integration should be 
performed over xs- If Cx does not depend on xs, in the final expression Eq. (|T8|) Sioc{crx/lu) 
appears as an additional factor at the oscillatory term. For ^ 1^ one has Sioc ~ 1 and 
Eq. (1481) reproduces the original result. Another way to take into account the localization 
of pions is to perfom the averaging of the oscillatory term in Eq. (ITS]) over the interval cr^. 

In the case cr^ ~ Idecay one expects some interplay between exponential decay and 
decoherence due to finite cXa; jS]. In this case in the limit ax 3> Idecay one should obtain again 
our original result. 




(47) 



Pfifi{Eu, L — x) 



P + 5';oc(o-x'Ai/)- sin^ 26'24 cos[^r(L - Xs)], 



(48) 
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4 MINOS data: Bounds and hint 



The number of the NC events with a given total energy of hadrons, E^, equals 

NiE,) = I dE^dxF^iE^)'^^^^^, (49) 

where cr{E^, Eh) is the cross-section of production of hadrons with energy E^ by neutrinos 
of energy E,^, y = E^/Ey is the inelasticity and x is the Bjorken variable. In what follows 
we will identify Eh with the reconstructed energy Ereco of [IJ. We assume that a{Ey,Eh) 
is well approximated by the deep inelastic scattering cross-section of neutrinos off nucleons. 
We use the MSTW 2008 set for the parton distribution functions and the reported spectrum 
for the NuMI beam (see [TU] for instance). 
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Figure 2: The ratio rTVD of events at the ND with and without oscillations as a function 
of the deposited energy for Am|2 = 0.5 eV^(circles), 1 eV^(squares), 2 eV^ (rhombuses) , 4 
eV^(upside triangles), 8 eV^(downside triangles). We take sin^ 26^24 = 0.1. 

Let N^{Eh) be the theoretical number of events computed according to Eq. fj49|) for 
Fy = Fy, i.e., as if there were no oscillations. We now define the ratio of events in the ND, 
with and without oscillations: 

N{Eh) 



N\Eh) ■ 



(50) 
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Figure 3: Predictions for the number of events at the ND (left) and FD 
oscillations. Solid histograms show predictions without 



right) with and 
- oscillations, 



The other histograms correspond to the case of z/^ — z/^ oscillations with 



without z/^ — z/< 
N% and iVO^ 

sin^ 26^24 = 0.15 and different values of Aml2: 0.5 eV^ (dotted), 2 eV^ (short dashed), 4 eV^ 
(long dashed), see Eqs. f lSTj) and fl53|) . The dots with the error bars represent the MINOS 
data points. 



The ratio tatd as a function of the hadron ("reconstructed") energy is shown in Fig. |2j Its 
behaviour with varying Am^g can be separated in three regions: 

1. For Am42 < 2 eV^ the peak of the spectrum of the incident neutrinos, E^, ~ 3.3 GeV, 
is in the region of small oscillation effect, (see Fig. [1]), and tat/) ~ 1 for most of the Eh 
range. The oscillation minimum is below Ey 2 GeV and, therefore, there is some 
depletion of events with energy E^ < 2 GeV. Nevertheless, r^n remains bigger than 
P for the whole energy range. 

2. For 2 eV^ < Am42 < 8 eV^ the dip in the oscillation probability falls in the region of 
energies for which the neutrino beam has maximum intensity. There is a considerable 
depletion of events for neutrino energies below 5 GeV and t^td < -P in this region. 
We will see below that, in this case, an excess of events should be observed in the FD 
compared to MINOS simulations. 

3. For Am|2 > 8 eV^ the first oscillation minimum is above the peak of the neutrino 
spectrum and the averaging is substantial for low Eh- Still, for Am|2 < 15 eV^ there 
are noticeable deviation from P below Eh = 5 GeV as MINOS becomes sensitive to 
higher oscillation minima. 

The Monte Carlo prediction N^jj for the number of events in the ND without oscilla- 
tions, (see Fig. 1 in [1]), needs to be corrected when the oscillations into sterile neutrino 
are present as 

NnD = VND N'^ND ■ (51) 
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At 

Am^^ = (1 - 3) eV^ sin^ 2^24 = 0.15 - 0.20 , (52) 

the oscillations can explain the deficit of events in the ND in the low energy bins (see 
Fig. [3]). This deficit is not yet statistically significant and for ND it is below the systematic 
errors. However the oscillations at ND must be taken into account when oscillation effect 
"propagates" from ND to FD. 

In Fig. [3] we show distributions of events over the reconstructed hadron energy Eh in 
the ND for different values of Aml2 (Eq. (53)). Also shown is the distribution without 
oscillations N^jj from [T] and experimental points. 

The MINOS collaboration predicts the number of events in the FD, Npjj, extrapolating 
the experimentally measured spectrum at the ND. In pQ it is assumed that oscillation 
effects are absent in the ND but the effect of the usual oscillations at the FD is included 
(see Fig. 2 in p]). The averaged z/^ — Ug oscillations in the FD modify this prediction: 
Npj^ — > PNpj-). However, when oscillations at the ND are taken into account, the latter 
must be corrected consistently with Eq. fl5T]) : 

Nfd = ^ N^fd- (53) 

Thus, the factor 

R{Eh) = (54) 
rND 

is the conversion factor for the flux predicted at the FD with and without of i/^ — z/^ oscil- 
lations. 

The behavior of R{Eh) can be straightforwardly inferred from Fig. [2j Again, there are 
three different cases: 

1. Averaged effect in ND: R{Eh) ~ 1. This happens at low energies for large Aml2 
when one has the averaged oscillation effect in the ND. The MINOS prediction for the 
number of events at the FD is unchanged by — oscillations. 

2. Strong suppression in ND: R{Eh) > 1. The neutrino energies in the peak of the 
neutrino flux are around the flrst oscillation minimum for 2 eV^ < Amj^ < 8 eV^. In 
this case one predicts some excess of events in the FD in comparison to the MINOS 
extrapolation which can partially explain the observed excess (see Fig. E]) 

3. Weak effect in the ND: R{Eh) < 1. This regime is realized for neutrino energies above 
the flrst oscillation minimum for all mass squared differences. It corresponds to small 
or no oscillation effect at the ND. 
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In Fig E] we show the distribution of events in the reconstructed hadron energy, Eh 
in the FD for different values of Aml2 (Eq. (55)) Also shown is the distribution without 
oscillations iV^^ from [Ij and experimental points. 

We consider next how the bounds on the — mixing are modified for different 
Am42. For Am\2 < 0.5 eV^ the decoherence is negligible, the oscillation effects in ND are 
significant only at low energies and therefore the MINOS limit in Eq. ([1]) is approximately 
valid. For Aml2 > 0.5 eV^ the limit should be modified. The strongest modification is for 
Am42 = (2 — 8) eV^, when the first oscillation minimum is at neutrino energies E^, that 
correspond to the peak in the spectrum of events. Notice that Eq. flS^ corresponds to 
^24 ~ (11 ^ 13)° which is substantially larger than the limit in Eq. ([T]). 

We stress that the calculations performed in this section are meant for illustration 
purposes. For precise quantitative results one should perform a complete MC simulation of 
events without the simplifications we made. 

5 Conclusions 

In this work we considered the oscillation effects in the ND of the MINOS experiment. The 
MINOS setup realizes an interesting situation of partial decoherence of the neutrino state 
at the production, when the energy splitting of two mass eigenstates is comparable with the 
energy uncertainty of the initial state (the width of pion). Decoherence leads to suppression 
of the oscillation effect at high energies (above the oscillation minimum), to the shift of 
oscillation minimum (dip) to low energies and to suppression of the depth of oscillations. 
The suppression becomes stronger at low energies. In general, the effect of decoherence 
should be taken into account for all experiments that perform searches for sterile neutrinos 
with 1 eV mass using neutrino beams from pion decays. 

The MINOS bounds, Eq. ([1]) remain valid for Am^2 < 0-5 eV^ For Am^2 > 0-5 eV^ the 
oscillation effect in the ND should be taken into account. It should noticeably weaken the 
bounds for Am|2 > 2 eV^. An estimation of this effect has been performed in this paper; a 
full quantitative analysis should be done by the MINOS collaboration. For Am42 > 15 eV^ 
the bound disappears and MINOS is insensitive to the z/^ — z/^ oscillations. 

The MINOS data might actually provide a hint of oscillations with Am\2 ~ (1 — 3) 
eV^. The oscillations can explain some deficit of signal in the low energy bins of the ND as 
compared to the Monte Carlo prediction, and the excess of events in the FD. 
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Appendix. Relation between the neutrino and pion en- 
ergies 



Using the Lorentz transformation to the laboratory frame we have for neutrino momenta: 

Py = po sin 6^0, Px ^ Po77r(cos6'o + (55) 

where 60 is the neutrino angle of emission in the pion rest frame; 7^ = ET^/m.,^ is the Lorentz 
factor of the pion and is the pion velocity. The neutrino energy in the laboratory frame, 
= -x/Px + Vyi ill the aproximation of ^ 1 can then be written as 

Ey ^ ^o77r(cos 00 + (3^) ^ £^o77r(cos 610 + 1). (56) 

Therefore, 

a^— ((coseo) + l)- (57) 

To evaluate the effective value of the angle (cos^o) let us consider the real experimental 
setup. According to Eq. (l55i) . the angle between the neutrino and pion in the laboratory 
frame equals 

sin^^o 
7^(cos6'o + 1) 

Consequently 

1 — 6 

cos 00 = b = T^tan^e^. (59) 

The maximal value of tan 6^, is given by 

tane'""^ = ^— (60) 
L — X 

with r 1 m being the radius of the fiducial zone of the ND and L is the distance from 
the target to the ND. It equals 2.5 x 10~^ for x = Ip and and 10"'^ for x = L. So, for the 
typical pion energy E^^ = 10 GeV we find the effective distance x = 300 m and b = 0.01, 
and therefore cos 6*0 ~ 0.98. We have checked that results do not change substantially for 
different values of a. 
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